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ji^ ■ Abstract. The double polarization (beam-recoil) observables Ox and Oz have been measured for the reac- 

^ ' tion 7p -^ K^ A from threshold production to E^ ~ 1500 MeV. The data were obtained with the linearly 

polarized beam of the GRAAL facility. Values for the target asymmetry T could also be extracted despite 
the use of an unpolarized target. Analyses of our results by two isobar models tend to confirm the necessity 
to include new or poorly known resonances in the 1900 MeV mass region. 

PACS. 13.60.Le Meson production - 13.88.-fe Polarization in interactions and scattering - 25.20.Lj Pho- 
toproduction reactions 

1 Introduction improvement by measuring cross sections with unprece- 
dented precision for a large number of channels but they 

A detailed and precise knowledge of the nucleon spec- also allowed a qualitative leap by providing for the first 

troscopy is undoubtedly one of the cornerstones for our un- time high quaUty data on polarization observables. It is 

derstanding of the strong interaction in the non-perturbative well known - and now well estabHshed - that these vari- 

regime. Today's privileged way to get information on the ables, being interference terms of various multipoles, bring 

excited states of the nucleon is light meson photo- and unique and crucial constraints for partial wave analysis, 

electroproduction. The corresponding database has con- hence facilitating the identification of resonant contribu- 

siderably expanded over the last years thanks to a com- tions and making parameter extraction more reliable, 

bined effort of a few dedicated faciUties worldwide. Not ^ , ■ ■ t^+ x , , • rr 

only did the recent experiments brought a quantitative ^rom this perspective, K+A photoproduction offers 

unique opportunities. Because the A is a self-analyzing 

Send offprint requests to: Ileres@lpsc.in2p3.fr particle, several polarization observables can be "easily" 
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measured via the analysis of its decay products. As a 
consequence, this reaction aheady possesses the richest 
database with results on the differential cross section Jj- 
[4], two single polarization observables (S and P) [2]-[6] 
and two double polarization observables {Cx and Cz) re- 
cently measured by the CLAS collaboration [7]. On the 
partial wave analysis side, the situation is particularly en- 
couraging with most models concluding to the necessity of 
incorporating new or poorly known resonances to repro- 
duce the full set of data. Some discrepancies do remain 
nonetheless either on the number of used resonances or 
on their identification. To lift the remaining ambiguities, 
new polarization obervables are needed calling for new ex- 
periments. 

In the present work, we report on first measurements 
of the beam-recoil observables Ox and Oz for the reaction 
-fp —i K^A over large energy (from threshold to 1500 
MeV) and angular (0^^ = 30 - 140°) ranges. The target 
asymmetry T, indirectly extracted from the data, is also 
presented. 



2 Experimental set-up 

The experiment was carried-out with the GRAAL facility 
(see [8] for a detailed description), installed at the Euro- 
pean Synchrotron Radiation Facility (ESRF) in Grenoble 
(France). The tagged and linearly polarized 7-ray beam is 
produced by Compton scattering of laser photons off the 
6.03 GeV electrons circulating in the storage ring. 

In the present experiment, we have used a set of UV 
lines at 333, 351 and 364 nm produced by an Ar laser, 
giving 1.40, 1.47 and 1.53 GeV 7-ray maximum energies, 
respectively. Some data were also taken with the green 
line at 514 nm (maximum energy of 1.1 GeV). 

The photon energy is provided by an internal tagging 
system. The position of the scattered electron is measured 
by a silicon microstrip detector (128 strips with a pitch of 
300 /im and 1000 /im thick). The measured energy resolu- 
tion of 16 MeV is dominated by the energy dispersion of 
the electron beam (14 MeV - all resolutions are given as 
FWHM). The energy calibration is extracted run by run 
from the fit of the Compton edge position with a precision 
of -lO/zm, equivalent to AE^/E^ ~ 2 x lO""* (0.3 MeV 
at 1.5 GeV). A set of plastic scintillators used for time 
measurements is placed behind the microstrip detector. 
Thanks to a specially designed electronic module which 
synchronizes the detector signal with the RF of the ma- 
chine, the resulting time resolution is «100 ps. The coin- 
cidence between detector signal and RF is used as a start 
for all Time-of-Flight (ToF) measurements and is part of 
the trigger of the experiment. 

The energy dependence of the 7-ray beam polarization 
was determined from the Klein-Nishina formula taking 
into account the laser and electron beam emittances. The 
UV beam polarization is close to 100% at the maximum 
energy and decreases smoothly with energy to around 60% 
at the KA threshold (911 MeV). Based on detailed studies 
[8], it was found that the only significant source of error 




Fig. 1. Schematic view of the LA7RANGE detector: BGO 
calorimeter (1), plastic scintillator barrel (2), cylindrical MW- 
PCs (3), target (4), plane MWPCs (5), double plastic scintil- 
lator hodoscope (6) (the drawing is not to scale). 



for the 7-ray polarization P^ comes from the laser beam 
polarization [5P^/P^=2%). 

A thin monitor is used to measure the beam flux (typ- 
ically 10*^ 7/s). The monitor efficiency (2.68±0.03%) was 
estimated by comparison with the response at low rate of 
a lead/scintillating fiber calorimeter. 

The target cell consists of an aluminum hollow cylin- 
der of 4 cm in diameter closed by thin mylar windows 
(100 /im) at both ends. Two different target lengths (6 and 
12 cm) were used for the present experiment. The target 
was filled by hquid hydrogen at 18 K (p « 7 10~^ g/cvc?). 

The 47r LA7RANGE detector of the GRAAL set-up 
allows to detect both neutral and charged particles (fig. 
[ij. The apparatus is composed of two main parts: a central 
one (25° < 6* < 155") and a forward one {9 < 25°). 

The charged particle tracks are measured by a set of 
MultiWire Proportional Chambers (MWPC) (see [5] for a 
detailed description). To cover forward angles, two plane 
chambers, each composed of two planes of wires, are used. 
The detection efficiency of a track is about 95% and the 
average polar and azimuthal resolutions are 1.5° and 2°, 
respectively. The central region is covered by two coaxial 
cylindrical chambers. Single track efficiencies have been 
extracted for 7r°p and n^n reactions and were found to be 
>90%, in agreement with the simulation. Since this paper 
deals with polarization observables, no special study was 
done to assess the efficiency of multi track events. Angular 
resolutions were also estimated via simulation, giving 3.5° 
in 9 and 4.5° in (p. 

Charged particle identification in the central region is 
obtained by dE/dx technique thanks to a plastic scintil- 
lator barrel (32 bars, 5 mm thick, 43 cm long) with an 
energy resolution «20%. For the charged particles emit- 
ted in the forward direction, a Time-of-Flight measure- 
ment is provided by a double plastic scintillator hodoscope 
(300x300x3 cm'^) placed at a distance of 3 m from the 
target and having a resolution of «600 ps. This detector 
provides also a measure of the energy loss dE/dx. Energy 
calibrations were extracted from the analysis of the 7r°p 
photoproduction reaction while the ToF calibration of the 
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forward wall was obtained from fast electrons produced in 
the target. 

Photons are detected in a BGO calorimeter made of 
480 {159x32ip) crystals, each of 21 radiation lengths. They 
are identified as clusters of adjacent crystals (3 on aver- 
age for an energy threshold of 10 MeV per crystal) with no 
associated hit in the barrel. The measured energy resolu- 
tion is 3% on average (^^=200-1200 MeV). The angular 
resolution is 6° and 7° for polar and azimuthal angles, 
respectively {E-^ > 200 MeV and ltarget=i cm). 



3 Data analysis 

3.1 Channel selection 

For the present results, the charged decay of the A {A ^ 
pn^ , BR=63.9%) was considered and the same selection 
method used in our previous publication on KA photo- 
production [5] was applied. Only the main points will be 
recalled in the following. 

Only events with three tracks and no neutral cluster 
detected in the BGO calorimeter were retained. In the 
absence of a direct measurement of energy and/or mo- 
mentum of the charged particles, the measured angles {9, 
If) of the three tracks were combined with kinematical 
constraints to calculate momenta. Particle identification 
was then obtained from the association of the calculated 
momenta with dE/dx and/or ToF measurements. 

The main source of background is the jp —> ptt+tt^ 
reaction, a channel with a similar final state and a cross 
section hundred times larger. Selection of the KA final 
state was achieved by applying narrow cuts on the follow- 
ing set of experimental quantities: 

. Energy balance. 

. Effective masses of the three particles extracted from 
the combination of measured dE/dx and ToF (only at 
forward angles) with calculated momenta. 



a 



. Missing mass m^p_ji+ evaluated from E^, 6k (mea- 
sured) and pk (calculated). 

For each of these variables, the width a of the corre- 
sponding distribution (Gaussian-like shape) were extracted 
from a Monte-Carlo simulation of the apparatus response 
based on the GEANT3 package of the CERN library 

To check the quality of the event selection, the dis- 
tribution of the A decay length was used due to its high 
sensitivity to background contamination. 

Event by event, track information and A momentum 
were combined to obtain the distance d between the re- 
action and decay vertices. The A decay length was then 
calculated with the usual formula ciyi = d/ {(3a * 1a)- Fig. 
[21 shows the resulting distributions for events selected with 
all cuts at ±2(7 (closed circles) compared with events with- 
out cuts (open circles). These spectra were corrected for 
detection efficiency losses estimated from the Monte-Carlo 
simulation (significant only for ct>5 cm). It should be 
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Fig. 2. Reconstructed A decay length spectrum after all se- 
lection cuts (closed circles) for events with at least two tracks 
in the cylindrical chambers. The solid line represents the fit 
with an exponential function a * exp{—ct/cT) where a and cr 
are free parameters. The second distribution (open circles) was 
obtained without applying selection cuts. It corresponds to the 
main background reaction {■yp —^ pn'^Tv') which, as expected, 
contributes only to small ct values. 



noted that the deficit in the first bins is attributed to fi- 
nite resolution effects not fully taken into account in the 
simulation. 

The first spectrum was fitted for ct>l cm by an ex- 
ponential function a * exp{—ct/cT) with a and ct as free 
parameters. The fitted ct value (8.17±0.31 cm) is in good 
agreement with the PDG expectation for the A mean free 
path (cr^=7.89cm) [9]. 

By contrast, the spectrum without cuts is dominated 
by pTT'^TT" background events. As expected, they contribute 
mostly to small ct values (<2-3 cm), making the shape of 
this distribution highly sensitive to background contami- 
nation. For instance, a pronounced peak already shows up 
when opening selection cuts at ±3cr. 

A remaining source of background, which cannot be 
seen in the ct plot presented above, originates from the 
contamination by the reaction 7p -^ K~^ S'^ . Indeed, events 
where the decay photon is not detected are retained by 
the first selection step. Since these events are kinemati- 
cally analyzed as KA ones, most of them are nevertheless 
rejected by the selection cuts. From the simulation, this 
contamination was found to be of the order of 2%. 

As a further check of the quality of the data sample, 
the missing mass spectrum was calculated. One should 
remcnrbcr that the missing mass is not directly measured 
and is not used as a criterion for the channel identifica- 
tion. The spectrum presented in fig. [3] (closed circles) is 
in fair agreement with the simulated distribution (solid 
line). Some slight discrepancies can nevertheless be seen 
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3.2.1 Formalism 
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Fig. 3. Distribution of the missing mass m^p_if+ recon- 
structed from measured E^, and 6k and calculated pk- Data 
after all selection cuts (closed circles) are compared to the sim- 
ulation (solid line). The expected contribution from the reac- 
tion 7p -^ K^ S° is also plotted (note that it is not centered 
on the S'^ mass due to kinematical constraints in the event 
analysis). The vertical arrow indicates the A mass. 



in the high energy tail of the spectra. The simulated miss- 
ing mass distribution of the contamination from the 7p -^ 
K^E^ reaction, also displayed in fig. [Sj clearly indicates 
that such a background cannot account for the observed 
differences. Rather, these are attributed to the summa- 
tion of a large number of data taking periods with vari- 
ous experimental configurations (target length, wire cham- 
bers, green vs UV laser line, ...). Although these config- 
urations were implemented in corresponding simulations, 
small imperfections (misalignments in particular) could 
not be taken into account. 

To summarize, thanks to these experimental checks, 
we are confident that the level of background in our se- 
lected sample is limited. This is corroborated by the sim- 
ulation from which the estimated background contamina- 
tion (multi-pions and K'^ S'^ contributions) never exceeds 
5% whatever the incident photon energy or the meson re- 
coil angle. 



For a linearly polarized beam and an unpolarized target, 
the differential cross section can be expressed in terms of 
the single polarization observables S, P, T (beam asym- 
metry, recoil polarization, target asymmetry, respectively) 
and of the double polarization observables Ox , Oz (beam- 
recoil), as follows [To]: 



[1 — P-yS cos2(p^ 



da 


1 (da 


PUn- 


2\dn 


+ax'P-yOx sin 2ip^ 




ay, {P - P^T cos 2(p^) 




-\-az'P^Oz sin2(/3-y] 





(1) 



Pf is the density matrix for the lambda final state and 
{da/df2)o the unpolarized differential cross section. The 
Pauli matrices a^'y.z' refer to the lambda quantization 
axes defined by z' along the lambda momentum in the 
center-of-mass frame and y' perpendicular to the reaction 
plane (fig. S]) . Py is the degree of linear polarization of the 
beam along an axis defined by n = a; cos ip^ +y sin ip^ ; the 
photon quantization axes are defined by z along the pro- 
ton center-of-mass momentum and y=y' (fig.H]). We have 
'Pf = flab — f, where fiab and (/? are the azimuthal angles 
of the photon polarization vector and of the reaction plane 
in the laboratory axes, respectively (fig.[S]). 

The beam-recoil observables Cx and Cz measured by 
the CLAS collaboration with a circularly polarized beam 
[7] were obtained using another coordinate system for de- 
scribing the hyperon polarization, the z' axis being along 
the incident beam direction instead of the momentum of 
one of the recoiling particles (see fig. H]). Such a non- 
standard coordinate system was chosen to give the re- 
sults their simplest interpretation in terms of polarization 
transfer but implied the model calculations to be adapted. 
To check the consistency of our results with the CLAS val- 
ues (see sect. 14. 1|) . our Ox and Oz values were converted 
using the following rotation matrix: 



O'^ = -Ox cos 9 cm - Oz sin 9 cm 
01 = Ox sin 6 cm ~ Oz cos 9cm 



(2) 



It should be noted that our definition for Ox and Oz 
(eq. [T]) has opposite sign with respect to the definition 
given in the article [llj , which is used in several hadronic 
models. We chose the same sign convention than the CLAS 
collaboration. 

For an outgoing lambda with an arbitrary quantization 
axis n' the differential cross section becomes: 



3.2 Measurement of Or, Oz and T 



As will be shown below, the beam-recoil observables Ox 
and Oz, as well as the target asymmetry T, can be ex- 
tracted from the angular distribution of the A decay pro- 
ton. 



dil 



da 



(3) 



where Pa is the polarization vector of the lambda. If the 
polarization is not observed, the expression for the differ- 
ential cross section reduces to: 



da 

dn 



Tr 



Pf 



da 
dn 



(4) 
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da 

dn 



{cpiab = 90°) 



da 

dn 



[l + P^Scoa2ip] (7) 



The beam asymmetry values S published in [S] were ex- 
tracted from the fit of the azimuthal distributions of the 
ratio: 



Njifilab = 90°) - Njifiiab = 0°) 
Ni^lab = 90°) + N{^lab = 0°) 



= P^Scos2ip (8) 



3.2.2 A polarization and spin observables 



Fig. 4. Definition of the coordinate systems and polar angles 
in the center-of-mass frame (viewed in the reaction plane) . The 
[x' ,y' ,z'] system is used to specify the polarization of the out- 
going A baryon: z' is along the A momentum and y' perpendic- 
ular to the reaction plane. The [a;,y,S] system is used to specify 
the incident photon polarization: z is along the incoming pro- 
ton momentum and y identical to y' . The polar angle 9cm of 
the outgoing K^ meson is defined with respect to the incident 
beam direction ziat- [x'c,y'c,z'c] is the coordinate system cho- 
sen by the CLAS collaboration for the A polarization. The x'^ 
and z'c axes are obtained from x' and z' by a rotation of angle 

71" -I- 6cm- 




Fig. 5. Definition of the coordinate systems and azimuthal an- 
gles in the center-of-mass frame (viewed perpendicularly to the 
beam direction). The [xiab,yiab,ziab] system corresponds to the 
laboratory axes with ziab along the incident beam direction. 
The [x,y,z] system, used to define the incident photon polar- 
ization, has its axes x and y along and perpendicular to the 
reaction plane (azimuthal angle ip), respectively. The polariza- 
tion of the beam is along n (azimuthal angle ^piab)- The two 
beam polarization states correspond to tpiab ~ O" (horizontal) 
and ipiab = 90° (vertical) [ipiab = 1/37 + V^)- 



which leads to: 



da 



da 

dn 



[1- P^Ucos2(p^] 



(5) 



For horizontal {ipiab = 0°) and vertical {(piab = 90°) pho- 
ton polarizations, the corresponding azimuthal distribu- 
tions of the reaction plane arc therefore: 



l^'«-»°>^(S)['-^^^™^''i 



(6) 



The components of the lambda polarization vector de- 
duced from eqs.[T]to[5]are: 



px ,z 



pV 

^A 



1 — P^S coa2(p~, 



P - P^T cos2(p^ 

1 — PjS cos2(pj 



(9) 



(10) 



These equations provide the connection between the A 
polarisation Pa and the spin observables S, P, T, Ox 
and Oz. 

Integration of the polarization components over the 
azimuthal angle ip of the reaction plane writes: 



<P\> 






(11) 



where i stands for x', y' or z' . 

When integrating over the full angular domain, the 
averaged x' and z' components of the polarization vec- 
tor vanish while the y' component is equal to P. On the 
other hand, when integrating over appropriatly chosen an- 
gular sectors, all three averaged components can remain 
different from zero. For horizontal and vertical beam po- 
larizations, the following expressions are obtained when 
considering the four particular ip domains defined here- 
after 12j (recalling ip^ = ipiab — 'p)'- 

. Si ^ [7r/4,37r/4] U [57r/4, 77r/4]: 

< pi > {ipiab - 0°) = [P^ + 2P^T)/(7r + 2P^S) 

< pi > {ipiab - 90°) = [Pt: - 2P-,T)/{'K - 2P^S) 

. 5-2 = [-7r/4,7r/4] U [37r/4, 57r/4]: 

< pi > {ifiab = 0°) = (Ptt - 2P^T)/(7r - 2P^i:) 

< pi > (ipiab - 90°) = (Ptt + 2P^T)/(^ + 2P^i:) 

. 53 = [0,7r/2]U[^,37r/2]: 

< Pi''' > {p>lab = 0°) = -2P^Ox^,/tt 

< Pi''' > {ipiab = 90°) = +2P^Ox,z/7T 



S4 = [tt/2, tt] U [37r/2, 2tt] 



<P1 

<p: 



> {flab 

> {flab 



0°) = 
90°) 



^2P^Ox,z/tt 
^2P^Ox,zh 
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It should be noted that these four sectors cover the full ^p 
range. In the following, these different combinations of ^p 
sectors and polarization states will be labelled by the sign 
plus or minus appearing in the corresponding expressions 



for < P\ >. 



3.2.3 Decay angular distribution 



3.2.4 Experimental extraction 

As for U and P, the observables Ox, Oz and T were ex- 
tracted from ratios of the angular distributions, in order 
to get rid of most of the distorsions introduced by the 
experimental acceptance. 

Including the detection efficiencies, the yields mea- 
sured as a function of the proton angle with respect to 
the different axes write: 



In the lambda rest frame, the angular distribution of the 
decay proton is given by |13| : 

W{cosep)^^{l + a\PA\cosep) (12) 

where q;=0.642±0.013 [S] is the A decay parameter and 
6p the angle between the proton direction and the lambda 
polarization vector. 

From this expression, one can derived an angular dis- 
tribution for each component of Pa '■ 



1 



Wicos 01) = -{1 + aP\ cos el) 



(13) 



where 0!, is now the angle between the proton direction 
and the quantization axis i {x\ y' or z'). 

The components being determined in the A rest frame, 
a suitable transformation should be applied to calculate 
them in the center-of-mass frame. However, as the boost 
direction is along the lambda momentum, it can be shown 
that the polarization measured in the lambda rest frame 
remains unchanged in the center-of-mass frame 0- 

When integrating over all possible azimuthal angles (/?, 
the proton angular distribution with respect to the y'-axis 
simply writes: 



w{cos ey')^^{i + aP cos ey' 



(14) 



where P is the recoil polarization. Our P results pub- 
lished in [S] were determined directly from the measured 
up/down asymmetry: 



N{cos 0y' > 0) - iV(cos 0y' < 0) 1 
iV(cos6i^' > 0) +iV(cos6i^' < 0) ^ 2 



(15) 



When integrating over the different angular domains 
specified above (sectors 5*1 -I- 5*2 for y'-axis and 6*3 + S4 for 
x'-,z'-axes, appropriatly combined with the two beam po- 
larization), the proton angular distributions with respect 
to the three quantization axes can be written as follows: 



W±{cost 



') = ^(1±« 



IP O 



■COS 



qx ,Z 



(16) 



M^±(cos<) = i(l + a^^^^^^cos<) (17) 



TT ± 2P^E 



/ / 1 / / , , 9P O 

i^tX .Z ^ iwtX .Z / riX .Z \ / 1 i J^^X,Z ^x ,Z \ 

Nj_ ^ -^Nq^ e±{cos 9p' )[l±a ■ cos 9p' ) 



m 



iiV^±e±(cos0^')(l-Ha 



Pn±2R/r 

n ± 2P^S 



(18) 



cos 61^') (19) 



From the integration of the azimuthal distributions given 
by eqs. [6] and [7| over the different angular sectors, it can 
be shown that: 



(20) 



No+' 
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:iV-^ 
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TT 


+ 2P^S 


Ny'_ 


TT 


-2P^S 



(21) 



Assuming that the detection efficiencies do not depend on 
the considered ip sectors {e+{cos6p) — e_(cos0p) - the va- 
lidity of this assumption will be discussed later on) , we can 
then calculate the following sums and ratios from which 
the efficiency cancels out: 



Nl 



N'_ 



-{N^^^ +N^o-n^±i^os0;'n 



(22) 



Nl+Nl =- (iVo\ + Ny_ )e± (cos 0^' ) (1 + aP cos 6^' ) 



2N^''' 



Nl 



2Nl 



Nl 



= (1 



2P (1 

t ^COS^^ ' 



N\ 



y' 



N_ 



v' 



(1 



2PjZ: ., 



"^ Tr+2P^ S ^°^ "^P 



1 



aP cos 6p 



(23) 



(24) 



(25) 



To illustrate the extraction method of Ox, O^ and T, 
the Nj^ and N- experimental distributions together with 
their sums and ratios, summed over all photon energies 
and meson polar angles, are displayed in figs. [6] (x'-axis), 
[7] (z'-axis) and [5] (y'-axis). Thanks to the efficiency cor- 
rection given by the distributions iV+ + iV_ (figs. 1617181 -c). 
the ratios 2N+/{N+ + N_) (figs.EEIld), from which the 
efficiency drops out, exhibit the expected dependence in 
cos 9p and can be therefore fitted by the functions given in 
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the r.h.s. of eg s . [24l and [25l The known energy dependence 
of Py and the previously measured values for S and P [5] 
are then used to deduce Ox, Oz and T from the fitted 
slopes. 

The vahdity of the hypothesis e+(cos6'p) — e-{cos 6p) 
was studied via the Monte Carlo simulation in which a 
polarized A decay was included. The efficiencies e± calcu- 
lated from the simulation are presented in plots e) of figs. 
[6] to [8] and the ratios e-/e+ in plots f) (open circles). As 
one can see, for the y' case, this ratio remains very close to 
1 whatever the angle while, for x' and z' , the discrepancy 
from 1 is more pronounced and evolves with the angle. 
This shows that some corrections should be applied on 
the measured ratios 2N+/{N^ + N^) to take into account 
the non-negligible differences observed between e+ and e_. 
The correction factors, plotted in figs. I6|7l8l -f) (closed cir- 
cles), were calculated through the following expression: 



Cor = 



2NI 



/ Qe7i' V , 



2N1 



'■Ni + Nl'9'^"' "-m+Ni 



' sel 



(26) 



where gen and sel stand for generated and selected events. 
Since e± — {N±)sei/iN±)gen, it can be re-writtcn as: 



Cor^ 



2N\ 



2^m 



Nl 



■i— 



(27) 



The corrected distributions are displayed in the plots g) of 
figs. [5] to [51 After correction, as expected, the slope of the 
y' distribution is unaffected while the slopes of the x' and 
z' distributions are slightly modified. These distributions 



To illustrate this second extraction method, the cor- 
rected distributions, summed over all photon energies and 
meson polar angles, are displayed in figs. l6l7l -i) {x' , z'-axes) 
and[8]-h),i) (y'-axis). They were obtained by dividing the 
originally measured distributions (figs. l6l7l -h and[Hl-a,b) by 
the corresponding efficiency distributions (figs. I6I7U and 
[8]-e). In the y'-axis case, the two corrected spectra N±/e± 
were simultaneously fitted. 

This method gives results in good agreement with those 
extracted from the first method. Nevertheless, the result- 
ing x^ were found to be significantly larger (the global 
reduced-x^ values are given in figs. [H] to [5]- they are close 
to 1 for the first method and five to ten times larger for 
the second one). The first method, which relies upon ratios 
leading to an intrinsic first order efficiency correction, is 
less dependent on the simulation details and was therefore 
preferred. 

Three sources of systematic errors were taken into ac- 
count: the laser beam polarization {6P-y/Pj=2%), the A 
decay parameter a {da — 0.013) and the hadronic back- 
ground. The error due to the hadronic background was 
estimated from the variation of the extracted values when 
cuts were changed from ±2cr to ±2.5cr. Given the good 
agreement between the two extraction methods, no corre- 
sponding systematic error was considered. For the T ob- 
servable, the measured values for S and P being involved, 
their respective errors were included in the estimation of 
the uncertainty. All systematic and statistical errors have 
been summed quadratically. 



were again fitted to obtain the final values of O 
T. 

As the detection efficiencies and the correction factors 
calculated from the simulation depend on the input values 
of Ox, Oz and T, an iterative method was used. Three 
iterations were sufficient to reach stable values. 



Oz and 4 Results and discussions 



The complete set of beam-recoil polarization and target 

asymmetry data is displayed in figs. [9] to [151 These data 

cover the production threshold region (i?-y=911-1500 MeV) 

^ . , , , • -1 ,and a large angular range (6*?:"°" = 30 — 140"). Numeri- 

For a consistency check, an alternative extraction method ^^1 ^^1^^^ ^^^ ^.^^^^ -^ tables [2 to [1 Error bars are the 

was implemented.^ The angukr distributions were directly q^^dratic sum of statistical and systematic errors. 

corrected by the simulated efficiencies and fitted according 

to: 
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(28) 



(29) 



Ptt 



2P-.,T 



o+: 



2P^i:' 



2p^j: 



■M){m 



4.1 Observable combination and consistency check 

In pseudoscalar meson photoproduction, one can extract 
experimentally 16 different quantities: the unpolarized dif- 
ferential cross section {da/df2)o, 3 single polarization ob- 
servables (P, T, S), 4 beam-target polarizations {E, F, 
G, H), 4 beam-recoil polarizations (Cx, Cz, Ox, Oz) and 
4 target-recoil polarizations {Tx, Tz, Lx, Lz). The various 
spin observables are not independent but are constrained 
by non-linear identities and various inequalites [10] . [llj . 
[2], [TS]. In particular, of the seven single and beam-recoil 
polarization observables, only five are independent being 
related by the two equations: 



where A^"'" and e'"" stand for N{- cos Op) and e(- cos Op), 
respectively. This trick, used for the x' and z' cases, al- 
lows to combine the N+ and N^ distributions which have 
opposite slopes (eq. [T8|. 



cl + c1 + ol + ol = i 



rpl p2 __ y^2 



CzOx - CxOz ^T-PE 



(31) 



(32) 
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Fig. 6. Angular distributions for the decay proton in the 
lambda rest frame with respect to the a;'-axis: a) distribu- 
tion A''+; b) distribution A'^_; c) sum A'^+ + A'^_; d) ratio 
2N+/{N+ + A^_); e) efficiencies e+ (triangles) and e_ (circles) 
calculated from the simulation; f) ratio e-/e+ (open circles) 
and correction factor Cor (closed circles) given by ea. l26l calcu- 
lated from the simulation; g) ratio 2N+/{N+ + A'^_) corrected 
by the factor Cor; h) distribution iV+ + Afl"", with N'"" = 



Ni-cosOp); i) efficiency e+ + e™", with e'' 



e(— cosSp 



calculated from the simulation; j) distribution A^+ + N^^ cor- 
rected by the efficiency e+ -I- e!!"" . The solid line in d) and g) 
represents the fit by the (linear) function given in the r.h.s. 
of eg. 1241 The solid line in j) represents the fit by the (linear) 
function given in the r.h.s. of eq. 1281 The reduced-x^ and the 
Ox value obtained from the fits are reported in d), g) and j). 

There are also a number of inequalities involving three of 
these observables: 



\T±p\ <i±z; 



p^ + ol + ol<i 



Oi 



Oi < 1 



(33) 
(34) 
(35) 



Fig. 7. Angular distributions for the decay proton in the 
lambda rest frame with respect to the z'-axis (all distributions 
as in fig. [Gjl . The reduced-x^ and the Oz value obtained from 
the fits are reported in d), g) and j). 



P^-hC: + Ct < 1 



r^ + c2 + c2 < 1 



(36) 



(37) 



These different identities and inequalities can be used 
to test the consistency of our present and previous mea- 
surements. They can also be used to check the compati- 
bility of our data with the results on Cx and Cz recently 
published by the CLAS collaboration 7 . 

Our measured values for S, P, T, Ox and Oz were 
combined to test the above inequalities. Equation 1311 was 
used to calculate the quantity C^ + Cl appearing in ex- 
pressions |36] and [37l The results for the two combinations 
T ± P| =F Z" of the three single polarizations are presented 
in fig. [121 The results for the quantities: 
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Fig. 8. Angular distributions for the decay proton in the 
lambda rest frame with respect to the j/'-axis: a) distribu- 
tion N+; b) distribution A''_; c) sum N+ + A''^; d) ratio 
2N+/{N+ +A''_); e) efficiencies e+ (triangles) and t- (circles) 
calculated from the simulation - they are symmetrical about 
0cm = 90° (we find edown /<^up — 1-03) ; f) ratio e-/e+ (open cir- 
cles) and correction factor Cor (closed circles) given by eq. 1261 
calculated from the simulation; g) ratio 2N+/{N+ -\- N-) cor- 
rected by the factor Cor; h) distribution A'^+ corrected by the 
efficiency e+; i) distribution N- corrected by the efficiency £_. 
The solid line in d) and g) represents the fit by the (non-linear) 
function given in the r.h.s. of eq. 1251 These distributions ex- 
hibit a linear behaviour since the overall recoil polarisation P 
is very low (the value extracted from the up/down asymmetry 
of the raw distribtion N+ -\- N- is -0.12). The solid line in h) 
and i) represents the simultaneous fit by the (linear) functions 
given in the r.h.s. of eas. l29l and l30l The reduced-x^ and the T 
value obtained from the fits are reported in d), g), h) and i). 



. {l+T^-P^-Ol- 02)l/2 ^ (^2 ^ ^2 ^ ^2)1/2^ 

. {1 + T^-E^-Ol- 01)'/^ = (P2 + cj + C2)i/2, 

which combine single and double polarization observables, 
are displayed in figs. fT3l and [T4l All these quantities should 
be < 1. The plotted uncertainties are given by the stan- 
dard error propagation. Whatever the photon energy or 
the meson polar angle, no violation of the expected in- 
equalities is observed, confirming the internal consistency 
of our set of data. 

Since all observables entering in eqs. [31] and [32] were 
measured either by GRAAL {S, P, T, Ox, Oz) or by 
CLAS (P, Cx, Cz - their P data were confirmed by our 
measurements [5]), the two sets of data can be there- 
fore compared and combined. Within the error bars, the 
agreement between the two sets of equal combinations 
{l+T^-E^-Ol-Oiy/^ (GRAAL) and (P^ + C^ + C^)!/^ 
(CLAS) is fair (fig. [T4)) and tends to confirm the previ- 
ously observed saturation to the value 1 of P = (P^ + 
Cx + C^y^"^, whatever the energy or angle. Fig. [15] dis- 
plays the values for the combined GRAAL-CLAS quan- 
tity CzOx ~ CxOz — P + PS. Within the uncertainties, 
the expected value (1) is obtained, confirming again the 
overall consistency of the GRAAL and CLAS data. 

It has been demonstrated [Tl] that the knowledge of 
the unpolarized cross section, the three single-spin observ- 
ables and at least four double-spin observables - provided 
they have not all the same type - is sufficient to determine 
uniquely the four complex reaction amplitudes. Therefore, 
only one additional double polarization observable mea- 
sured using a polarized target will suffice to extract un- 
ambiguously these amplitudes. 



4.2 Comparison to models 

We have compared our results with two models: the Ghent 
isobar RPR (Regge-plus-resonance) model [16]- [19] and 
the coupled-channel partial wave analysis developed by 
the Bonn-Gatchina collaboration 20j-[24 . In the follow- 
ing, these models will be refered as RPR and BG, respec- 
tively. The comparison is shown in figs. Olto llll 

The RPR model is an isobar model for KA photo- and 
electroproduction. In addition to the Born and kaonic con- 
tributions, it includes a Reggeized t-channel background 
which is fixed to high-energy data. The fitted database in- 
cludes differential cross section, beam asymmetry and re- 
coil polarization photoproduction results. The model vari- 
ant presented here contains, besides the known N* reso- 
nances (S'ii(1650), Pii(1710), Pi3(1720)), the Pi3(1900) 
state (** in the PDG [5]) and a missing P>i3(1900) reso- 
nance. This solution was found to provide the best overall 
agreement with the combined photo- and electroproduc- 
tion database. As one can see in figs. [9] to [11] the RPR 
prediction (dashed line) qualitatively reproduces all ob- 
served structures. Interestingly enough, the model best re- 
produces the data at high energy (1400-1500 MeV), where 
the Pi3(1900) and 1)13(1900) contributions are maximal. 

The BG model is a combined analysis of experiments 
with ttN , rjN, KA and KS final states. As compared 
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to the other models, this partial- wave analysis takes into 
account a much larger database which includes most of 
the available results (differential cross sections and po- 
larization observables). For the 7p -^ K^A reaction, the 
main resonant contributions come from the S'ii(1535), 
S'ii(1650),Pi3(1720),Pi3(1900) and Pii{1840) resonances. 
To achieve a good description of the recent Cx and Cz 
CLAS measurements, the ** Pi3(1900) had to be intro- 
duced. It should be noted that, at this stage of the anal- 
ysis, the contribution of the missing _Di3(1900) is signif- 
icantly reduced as compared to previous versions of the 
model. As shown in figs. [MTTj this last version (solid line) 
provides a good overall agreement. On the contrary, the 
solution without the Pi3(1900) (not shown) fails to repro- 
duce the data. 

More refined analyses with the RPR and BG models 
are in progress and will be published later on. Compari- 
son with the dynamical coupled-channel model of Saclay- 
Argonne-Pittsburgh [25]- [IJ has also started. 



5 Summary 

In this paper, we have presented new results for the re- 
action 7p —>■ K+A from threshold to E^ '- 1500 MeV. 
Measurements of the beam-recoil observables Ox , Oz and 
target asymmetries T were obtained over a wide angu- 
lar range. We have compared our results with two isobar 
models which are in reasonable agreement with the whole 
data set. They both confirm the necessity to introduce 
new or poorly known resonances in the 1900 MeV mass 
region (P13 and/or Di^). 

It should be underlined that from now on only one 
additional double polarization observable (beam-target or 
target-recoil) would be sufficient to extract the four helic- 
ity amplitudes of the reaction. 
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Fig. 9. Angular distributions of the beam recoil observable Ox- Data are compared with the predictions of the BG (solid line) 
and RPR (dashed line) models. 
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Fig. 10. Angular distributions of the beam recoil observable O^. Data are compared with the predictions of the BG (solid line) 
and RPR (dashed line) models. 
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Fig. 11. Angular distributions of the target asymmetry T. Data are compared with the predictions of the BG (solid line) and 
RPR (dashed line) models. 
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Fig. 12. Angular distributions of the quantities \T — P\ -{- E (closed circles) and |r + P| — X" (open circles). We should have 
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Fig. 13. Angular distributions of the quantities (P^ + O^ + O^)^/^ (stars), (r^ + O^ + O^)^/^ (circles) and (1 + T^ - P^ - 
Ox — Oz)^ = (-S'^ + C^ + C'^y (crosses). The first and third sets of data are horizontally shifted for visualization. All these 
quantities should be < 1 (inequalities 1341 135] and l37 |l . 
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Fig. 14. Angular distributions of the quantity (1 + T^ - T^ - O^ - O'iY^"^ = (P^ + Cl + Clf^^. This quantity should be < 1 
(inequality [36} . Comparison to the values {P^ + C^ + ClY published by the CLAS collaboration (open squares - energy in 
parentheses). Note that the O^ + O^ and C^ + C^ values are independent of the choice for the axes specifying the A polarization 
(see sect. l3.2Tl]) . 
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Fig. 15. Angular distributions of the quantity C^Ox — C^Oz ~ T + PS. This quantity is calculated using the Cx and Cz 
results published by the CLAS collaboration (energy in parentheses) combined with our Ox and Oz data converted by eq. [2] 
to have the same z axis convention and with our S, P and T measurements. The used CLAS data are those corresponding 
to the angles cos0cm=O.85, mean(0.65,0.45), mean(0.25,0.05), -0.15, mean(-0. 35,-0. 55) and -0.75. We should have the equality 
CzOx-CxOz -T + PS = (eq.[32l). 
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Table 1. Beam-recoil Ox values. 



ecm(°) 


S^=980 MeV 


ecm(°) 


E.-IQ27 McV 


ecAl 


E^-IQIA MeV 


Oc^n 


£^ = 1122 MeV 


31.3 


0.349 ± 0.150 


30.6 


0.425 ± 0.108 


31.2 


0.502 ± 0.103 


32.4 


0.570 ± 0.154 


59.1 


0.320 ± 0.255 


57.5 


0.408 ± 0.225 


57.5 


0.202 ± 0.140 


57.6 


0.179 ± 0.161 


80.7 


-0.094 ± 0.262 


81.7 


-0.085 ± 0.189 


81.0 


-0.190 ± 0.120 


80.6 


-0.365 ± 0.282 


99.8 


-0.464 ± 0.320 


99.8 


-0.477 ± 0.133 


99.7 


-0.552 ± 0.106 


99.2 


-0.522 ± 0.165 


118.9 


-0.490 ± 0.244 


119.0 


-0.723 ± 0.142 


119.3 


-0.621 ± 0.105 


119.9 


-0.795 ± 0.146 


138.6 


-1.028 ± 0.410 


139.5 


-0.304 ± 0.259 


138.8 


-0.163 ± 0.200 


139.3 


-0.341 ± 0.252 


e^r^n 


E--ini McV 


e^r^n 


S^=1222 MeV 


e^r^n 


S^=1272 MeV 


Ocran 


£^=1321 MeV 


34.1 


0.567 ± 0.122 


34.6 


0.294 ± 0.179 


35.8 


0.526 ± 0.130 


36.0 


0.599 ± 0.161 


58.9 


0.306 ± 0.148 


58.9 


0.310 ± 0.139 


59.2 


0.345 ± 0.140 


59.4 


0.304 ± 0.155 


80.5 


-0.319 ± 0.254 


80.4 


-0.193 ± 0.124 


80.6 


0.028 ± 0.155 


80.3 


-0.073 ± 0.142 


99.3 


-0.510 ± 0.175 


98.9 


-0.548 ± 0.252 


99.2 


-0.232 ± 0.125 


99.2 


-0.383 ± 0.141 


119.4 


-0.347 ± 0.162 


119.1 


-0.615 ± 0.121 


119.9 


-0.395 ± 0.143 


119.9 


0.046 ± 0.114 


140.4 


-0.160 ± 0.234 


140.3 


0.116 ± 0.242 


140.8 


0.454 ± 0.187 


141.3 


0.323 ± 0.158 






Ocmn 


£^=1372 MeV 


Ocrnn 


£^=1421 McV 


dcm\) 


St- = 1466 MeV 






36.1 


0.552 ± 0.119 


35.7 


0.455 ± 0.144 


35.9 


0.307 ± 0.150 








59.5 


0.150 ± 0.130 


59.6 


-0.072 ± 0.160 


59.3 


0.172 ± 0.171 








80.1 


-0.168 ± 0.131 


80.3 


-0.303 ± 0.139 


80.0 


-0.270 ± 0.195 








99.4 


-0.276 ± 0.122 


99.7 


-0.190 ± 0.137 


99.7 


-0.096 ± 0.139 








120.4 


0.124 ± 0.138 


120.4 


0.076 ± 0.110 


120.8 


0.164 ± 0.145 








141.9 


0.636 ± 0.126 


142.8 


0.490 ± 0.205 


143.7 


0.905 ± 0.152 





Table 2. Beam-recoil Oz values. 



ec^n 


E^='dm MeV 


OcA") 


£;^=1027 MeV 


Ocr^n 


E^=IQ7A MeV 


Ocr^n 


£^=1122 MeV 


31.3 


0.581 ± 0.194 


30.6 


0.333 ± 0.110 


31.2 


0.285 ± 0.080 


32.4 


0.274 ± 0.124 


59.1 


0.956 ± 0.242 


57.5 


0.951 ± 0.216 


57.5 


0.674 ± 0.112 


57.6 


0.687 ± 0.127 


80.7 


0.754 ± 0.186 


81.7 


0.995 ± 0.154 


81.0 


1.003 ± 0.148 


80.6 


0.888 ± 0.244 


99.8 


1.139 ± 0.237 


99.8 


0.949 ± 0.140 


99.7 


1.140 ± 0.130 


99.2 


0.950 ± 0.144 


118.9 


0.841 ± 0.215 


119.0 


0.744 ± 0.162 


119.3 


0.996 ± 0.156 


119.9 


0.618 ± 0.197 


138.6 


-0.091 ± 0.597 


139.5 


-0.287 ± 0.415 


138.8 


0.427 ± 0.223 


139.3 


-0.162 ± 0.568 


Ocm\ ) 


£^=1171 McV 


6'cm(°) 


S^=1222 McV 


Ocr^O 


E^=1212 McV 


ecm(°) 


£^ = 1321 MeV 


34.1 


0.398 ± 0.093 


34.6 


0.291 ± 0.177 


35.8 


0.532 ± 0.087 


36.0 


0.554 ± 0.090 


58.9 


0.914 ± 0.128 


58.9 


0.678 ± 0.167 


59.2 


0.710 ± 0.119 


59.4 


0.904 ± 0.108 


80.5 


0.825 ± 0.123 


80.4 


0.485 ± 0.109 


80.6 


0.867 ± 0.112 


80.3 


0.767 ± 0.153 


99.3 


0.964 ± 0.175 


98.9 


1.025 ± 0.143 


99.2 


0.676 ± 0.188 


99.2 


0.734 ± 0.161 


119.4 


0.550 ± 0.190 


119.1 


0.426 ± 0.166 


119.9 


0.677 ± 0.166 


119.9 


0.409 ± 0.229 


140.4 


-0.055 ± 0.286 


140.3 


-0.162 ± 0.276 


140.8 


0.349 ± 0.272 


141.3 


-0.448 ± 0.217 






Ocmn 


£.^=1372 MeV 


OcA") 


£^=1421 MeV 


Oa^n 


St =1466 MeV 






36.1 


0.600 ± 0.084 


35.7 


0.384 ± 0.094 


35.9 


0.354 ± 0.095 








59.4 


0.784 ± 0.119 


59.6 


0.558 ± 0.185 


59.3 


0.814 ± 0.222 








80.1 


0.484 ± 0.112 


80.3 


0.322 ± 0.195 


80.0 


0.666 ± 0.332 








99.4 


0.419 ± 0.120 


99.7 


0.289 ± 0.134 


99.7 


-0.023 ± 0.192 








120.4 


0.019 ± 0.145 


120.4 


-0.313 ± 0.131 


120.8 


-0.432 ± 0.180 








141.9 


-0.072 ± 0.159 


142.8 


-0.085 ± 0.172 


143.7 


-0.461 ± 0.162 
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Table 3. Target asymmetry T values. 



6cm{°) 


S^=980 MeV 


ecm(°) 


E.-IQ27 MeV 


OcmO 


S^=1074 MeV 


ecmH 


£^=1122 MeV 


31.3 


-0.506 ± 0.156 


30.6 


-0.663 ± 0.112 


31.2 


-0.635 ± 0.096 


32.4 


-0.615 ± 0.118 


59.1 


-0.607 ± 0.206 


57.5 


-0.860 ± 0.190 


57.5 


-0.718 ± 0.120 


57.6 


-0.991 ± 0.237 


80.7 


-0.803 ± 0.185 


81.7 


-0.749 ± 0.139 


81.0 


-0.874 ± 0.141 


80.6 


-0.949 ± 0.239 


99.8 


-0.622 ± 0.166 


99.8 


-0.974 ± 0.121 


99.7 


-1.048 ± 0.140 


99.2 


-0.833 ± 0.153 


118.9 


-0.622 ± 0.187 


119.0 


-0.789 ± 0.133 


119.3 


-0.760 ± 0.102 


119.9 


-0.825 ± 0.165 


138.6 


-1.090 ± 0.341 


139.5 


-0.681 ± 0.359 


138.8 


-0.448 ± 0.203 


139.3 


-0.465 ± 0.296 


e^r^n 


£^=1171 McV 


e^r^n 


S^=1222 MeV 


OcraH 


S^=1272 MeV 


OcraH 


£^=1321 MeV 


34.1 


-0.715 ± 0.116 


34.6 


-0.858 ± 0.155 


35.8 


-0.773 ± 0.123 


36.0 


-1.064 ± 0.133 


58.9 


-0.869 ± 0.154 


58.9 


-0.874 ± 0.145 


59.2 


-0.827 ± 0.166 


59.4 


-0.910 ± 0.142 


80.5 


-0.850 ± 0.158 


80.4 


-0.871 ± 0.124 


80.6 


-0.979 ± 0.234 


80.3 


-0.716 ± 0.133 


99.3 


-0.659 ± 0.159 


98.9 


-0.690 ± 0.132 


99.2 


-0.707 ± 0.150 


99.2 


-0.576 ± 0.223 


119.4 


-0.669 ± 0.150 


119.1 


-0.675 ± 0.145 


119.9 


-0.125 ± 0.208 


119.9 


-0.281 ± 0.174 


140.4 


0.226 ± 0.316 


140.3 


-0.066 ± 0.249 


140.8 


0.482 ± 0.213 


141.3 


0.331 ± 0.198 






OcraH 


£^=1372 McV 


OcA") 


£^=1421 MeV 


Oa^n 


St =1466 MeV 






36.1 


-0.983 ± 0.104 


35.7 


-0.753 ± 0.112 


35.9 


-0.632 ± 0.127 








59.5 


-0.695 ± 0.113 


59.6 


-0.687 ± 0.159 


59.3 


-0.648 ± 0.166 








80.1 


-0.669 ± 0.123 


80.3 


-0.564 ± 0.131 


80.0 


-0.553 ± 0.185 








99.4 


-0.482 ± 0.175 


99.7 


-0.025 ± 0.157 


99.7 


0.190 ± 0.196 








120.4 


-0.104 ± 0.135 


120.4 


0.160 ± 0.112 


120.8 


0.785 ± 0.195 








141.9 


0.629 ± 0.147 


142.8 


0.859 ± 0.140 


143.7 


0.933 ± 0.175 





